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Abstract 

The input to the unrooted traveling repairman problem is an undirected metric graph and 
a subset of nodes, each of which has a time window of unit length. Given that a repairman 
can start at any location, the goal is to plan a route that visits as many nodes as possible 
during their respective time windows. A polynomial-time bicriteria approximation algorithm 
is presented for this problem, gaining an increased fraction of repairman visits for increased 
speedup of repairman motion. For speedup s, we find a 67/(5 -I- l)-approximation for s in the 
range 1 < s < 2 and a 47/s-approximation for s in the range 2 < s < 4, where 7 = 1 on 
£^ , tree-shaped networks and 7 = 2 + e on general metric graphs. 



1 Introduction 



When planning the route for an agent who needs to drop off supplies, make repairs, deliver impor- 
CN ! tant information, or other similar tasks, distance traveled may not be the only important metric. 

The usefulness of visiting a particular location may be heavily dependent on time. Consider the 
■ example of a client who schedules a window of time during which he or she will be home. An 

! appliance repairman will be able to perform repairs only during that time window. 

In the last decade, the algorithms community has achieved much progress with time-sensitive 
routing problems of this kind. These problems typically identify the locations to be visited and 
the cost of traveling between them as the nodes and edges, respectively, of a weighted, undirected 
graph. For example, the orienteering problem [U [21 IS] seeks to find a path that visits as 
^ ■ many nodes as possible before a global time deadline. The deadline traveling salesman problem [2] 

generalizes this problem further by allowing each location to have its own deadline. This problem 
can be generalized even further to the traveling repairman problem [21 (H [8j [T2l [1~9] by allowing 
each location to have a time window during which it must be visited to receive credit. 

In [12], we introduced the first polynomial-time algorithms that give constant approximations 
to the traveling repairman problem with unit-time windows whenever the underlying graph is a 
tree or a metric graph. In this version of the problem, a time window is identified with each service 
request whose location is given by the node with which it is associated. Let a repairman earn a 
specified profit for visiting a service request during its time window. A planned sequence of visits 
made by traveling at a given speed is called a service run. The goal of the traveling repairman 
problem is to plan a service run that maximizes profit. Unlike much of the preceding literature 
[21 H], we considered the unrooted version of the problem, in which the repairman may start at any 
time from any location and stop similarly. Although the unrooted problem is no harder than the 
rooted problem, which specifies a starting point, it is a fascinating and difficult problem in its own 
right. Both problems are NP-hard when the graph is a tree [12] and APX-hard for a general metric 
graph [13] . 
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As a counterpoint to the repairman problem, we also introduced the speeding deliveryman 
problem in [121 113] . with an alternative optimization paradigm, namely speedup. The input to 
the speeding deliveryman problem is the same as the input to the traveling repairman problem, 
but the goal is to find the minimum speed necessary to visit all service requests during their time 
windows and thus collect all profit. The deliveryman problem is similarly hard, and in [12\ [T3] we 
introduced polynomial-time approximation algorithms for it. 

In both the repairman and deliveryman problems, our algorithms from [121113] rely on trimming 
windows so that the resulting time windows are pairwise either identical or non-overlapping. To 
trim time windows, we first make divisions in time every .5 time units, starting at time 0. We call 
the time interval that starts at a particular division and continues up to the next division a period. 
Because we allow no window to start on a period boundary, each time window will completely 
overlap exactly one period and partially overlap the two neighboring periods. Trimming then 
removes those parts of each window that fall outside of the completely overlapped period. 

This process of trimming incurs a penalty in our approximations that is a reduction by a factor 
of 1/3 in the number of requests serviced by the repairman and an increase by a factor of 4 in 
the speedup needed by the deliveryman. Yet one might expect that a spectrum of performance is 
possible between these two extremes. If we have any particular speedup s greater than 1 but less 
than 4, we expect an increase in the number of serviced requests, proportionate in some sense to 
s. Indeed, this is so, as we show in this paper. 

The concept of speedup has been used in the scheduling community since its introduction in 
[15j . Subsequent work in [18] refined the technique and gave it the name resource augmentation 
as well as notation suitable for scheduling. Results in [3] have made notable strides by producing 
0(l)-speed 0(l)-approximation algorithms for some kinds of scheduling problems. In other words, 
by increasing the speed of the machines on which the jobs are run by a constant factor, the resulting 
schedule can achieve a measure of performance within a constant factor of the optimal schedule run 
on machines at a baseline speed. Although we do not use the same techniques or notation given 
in these papers, we note the connection between scheduling and time-constrained routing and the 
value of resource augmentation when finding approximation algorithms for each class of problem. 

For unit-time windows we present an algorithm that finds approximations parameterized by 
speedup s and a factor 7, where 7 = 1 for a tree and 7 is no more than 2 + e for a metric graph. 
These values for 7 are explained in Sect. [2] For s in the range 1 < s < 2, our algorithms find a 
67/(5 + l)-approximation. For s in the range 2 < s < 4, our algorithms find a 47/s-approximation. 
As long as unit speed is no greater than the speed needed for an optimal service run to collect 
all possible profit, our results will hold. Were the unit speed considerably faster than is needed 
to service all requests, the idea of offsetting the sub-optimality of an approximation by increasing 
speed beyond that of an optimal run becomes meaningless. 

Our analysis depends upon the interactions between speedup and trimming, but setting the 
starting point for partitioning periods at time as we did in [12J may not result in the most 
profitable run over trimmed windows. Especially in the case of non-integer speedups, careful study 
suggests that some particular starting point for periods may be better than others. If we repeatedly 
choose different starting points, run our repairman algorithm on trimmed windows, and always keep 
the best answer found so far, we can establish better bounds on the final approximation. Although 
there are an infinite number of points to start partitioning periods from, we can find a canonical 
set of starting points that are representative of all possible starting points, and the size of this set 
is linear in the number of service requests. 

Then, for each set of periods we use, we trim the windows and run a basic repairman algorithm. 
For each set of trimmed windows, our analysis depends on an ensemble of runs, each of which moves 
backward and forward along the path of an optimal run. The ensemble is composed of a number 
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of runs traveling at a speedup s faster than an optimal run travels. While it is too expensive to 
determine an actual optimal run unless P = NP, our analysis will rely only on a definition of each 
run in the ensemble in terms of the optimal run. One run in the ensemble might be effective in 
servicing requests serviced by an optimal run very early in their respective windows. Another run 
might be effective in servicing requests serviced by an optimal run very late in their respective 
windows. Together, these runs would service every request on trimmed windows that an optimal 
run would service on untrimmed windows. An optimal run on trimmed windows, then, would 
achieve at least the average service of the ensemble of runs taken as a whole. Likewise, the best 
run from among all possible choices of trimming would achieve at least the average service taken 
over all runs in the ensemble. 

Except for the factor of 7 that does not come from trimming, our algorithm finds a run which 
performs as well as this average. Our analysis of an ensemble guarantees that a certain fraction 
of profit is always collected by the run produced by our algorithm. Again, we know the ensemble 
of runs only symbolically, because we cannot efficiently know the instantiation of these runs unless 
P = NP. The challenge has been to devise ensembles so that our analysis will establish a good 
performance for our algorithm. 

Although we are the first to extend time window problems into the realm of speedup, much 
other related work is being done with time windows. For general metric spaces and general time 
windows together in the rooted problem, an 0(log 2 n)-approximation is given in [2]. An O(logL)- 
approximation is given in [6], for the case that all time window start and end times are integers, 
where L is the length of the longest time window. In constrast, a constant approximation is given 
in [8], but only when there are a constant number of different time windows. Following the initial 
publication of our work in [12J , an extension was given in [6] that gives an 0(log D)-approximation to 
the unrooted problem with general time windows, where D is the ratio of the length of largest time 
window to the length of the smallest. Poly logarithmic approximation algorithms to the directed 
traveling salesman problem with time windows have been given in [7] and [17j . The repairman 
problem with time windows has also been studied in the operations research community, as in |10j 
and where it is exhaustively solved to optimality. 

Note that the bicriteria approximation that we give will guarantee the repairman a greater 
fraction of the optimal profit by speeding up beyond the speed of the optimal service run. Although 
we are the first to identify techniques that allow for a trade-off between requests serviced and 
speedup, there are other bicriteria approximations for repairman problems. For example, a bicriteria 
approximation is presented in [2] that allows the repairman to service more requests if the time 
windows are all lengthened by some amount e. If the time window for request Sj starts at time t 
and ends at time t 1 , this lengthening would change each time window from [t, t') to [t, (1 + e)t'). 
In contrast, the effect of running at speedup s is equivalent to increasing all times by some factor, 
changing time window [t, t') to [st, st'). 

Before we explain our bicriteria algorithms in depth, we will summarize in Sect. [2] our results 
from [12j as well as recent improvements. In Sect. El we will explain the basic mechanisms needed 
for the averaging analysis used in our ensemble approach. In Sect. HI we will introduce the concept 
of canonical sets needed to make our algorithms practical for arbitrary values of speedup s. In 
Sect. 02 we will give the analysis that proves our approximation results for speedup s in the range 
2 < s < 4. Finally, in Sect. El we will give the more involved analysis proving approximation results 
for speedup s in the range I < s < 2. 
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2 Algorithms for Repairman without Speedup 



We continue to restrict our focus to the traveling repairman problem with unit-time windows. For 
this problem, our repairman approximation algorithms in [12t 113] trim the time windows so that 
the resulting windows can be partitioned into non-overlapping periods. These algorithms identify 
a variety of good paths within each separate period and then use dynamic programming to select 
and paste these paths together into a variety of longer good paths and ultimately a good service 
run for the whole problem. Optimal paths within each period can be found in polynomial time 
when the underlying graph is a tree, but only approximately optimal paths are found on a metric 
graph. 

To unify the notation between these approaches, we describe the approximation factor for finding 
paths within each period using 7, where 7 = 1 for a tree and 7 = 2 + e for a metric graph. Results 
from [12] clearly show that 7 < 5 for a metric graph, but this result has been improved to 2 + e 
because of results in [7J. To describe the running time for these repairman algorithms we bound 
the time by T(n), where T(n) is 0(n 4 ) for a tree and 0{n°( l l e )) for a metric graph, where the tree 
bound comes from results in [12] and the metric graph bound can be derived by augmenting those 
results with techniques from [7J. A full explanation of this improved result is available in |13j . 

3 The Ensemble Approach for Analyzing Performance 

The two results cited above give analysis for runs after trimming for the cases of no speedup (s = 1 
guaranteeing 1/(37) °f optimal profit) and full speedup (s = 4 guaranteeing I/7 of optimal profit 
achieved at the original speed). Our analysis cannot examine speeds lower than the reference speed 
of 1 (i.e, slowdown) because there is no guarantee that a slower run can ever earn any constant 
fraction of the profit of a faster run. Likewise, above a speedup of 4, the effects of trimming are 
completely offset. Thus, this range of speedups is the only range for which our bicriteria analysis 
seems practical. For this range 1 < s < 4, our analysis uses a number of different service runs that 
are defined by moving backward and forward along an optimal tour R*. We rely on averaging over 
a suitable ensemble of runs to establish that the run R generated by our algorithm does well. 

Let R* be an optimal service run at unit speed originally starting at time 0. In defining the 
service runs that we use in our analysis, we use the term racing to describe movement, forward and 
backward, along R* at a speedup of s times the speed at which R* travels. Please refer to Figure [1] 
for an example of runs racing at s = 2. For convenience, let 5 = l/(2s), which is the time needed 
to travel at speed s the distance traveled in one period at unit speed. 

One of the runs over which we average is service run A, defined as follows. Start run A at time 
at the location that R* has at time —.5. Then run A follows a pattern of racing forward along 
R* for 1 period, racing backward along R* for 1 — 25 periods, and then racing forward along R* 
for 25 periods. The positions reached by run A at the times 0, .5, 1, 1.5, and 2 are indicated in 
Figure [TJ Note that the pattern of movement for run A repeats every 2 periods. 

Define A R , the " reverse" of A, as follows. Start run A R at time at the location that R* has at 
time .5. Then run A R follows a repeating pattern of racing forward along R* for 25 periods, racing 
backward along R* for 1 — 25 periods, and then racing forward along R* for 1 period. Similarly, 
the positions reached by run A R at times 0, .5, 1, 1.5, and 2 are indicated in the same figure. 

If a service request p is serviced by a run R during the period that the time window of p has 
been trimmed into, we say that R covers p. We partition the set of requests whose windows are 
trimmed into a given period into three sets: Set T contains those requests serviced by R* in that 
same period, set E contains those requests serviced by R* in the preceding period, and set L 
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contains those requests serviced by R* in the following period. See Figure Q] for an example: Run A 
services requests in L for each period, and it services requests in T for every second period. Note, 
for example, that requests serviced on R* between times and .5 are serviced by A during that 
same period. 



Optimal 




-.5 



Figure 1: Runs A and A R with a speedup of 2 based on an optimal run. The pattern in A R perfectly 
complements A to give full coverage. Times are labeled on the optimal run as well as runs A and 
A R . Segments of each run are also designated L, T, and E depending on which subset of runs they 
add coverage to. The light vertical lines indicate shared positions on the runs, not identical times. 

We number the periods —1, 0, 1, 2, and so on by the integer multiples of .5 that give their 
starting times. The labels L, T, and E in Figure [T] make it clear that every period numbered 1 
and later has full coverage. The combination of runs A and A R ensures that the location of a 
given request is visited three times: during the same period, the previous period, and the following 
period that the request was visited by R*. Thus, we guarantee that our combination of runs visits 
the request no matter which period its window was trimmed into. However, it appears that periods 
— 1 and only have partial coverage. This apparent discrepancy is an endpoint anomaly. In fact, 
period —1 needs no coverage because no locations visited by R* are visited in this period. Likewise, 
there are no requests in set E serviced by R* during period because such requests would have 
been trimmed into period —1, which is impossible. 

Although there might be some superficial resemblance between the L, T, and E sets and the 
small margin and large margin cases discussed in [2], the approach in that paper iteratively runs 
specialized algorithms that progressively redefine the meaning of small margins and large margins. 
In contrast, the three sets that we use are fixed for any given trimming scheme, do not have separate 
algorithms tailored for each one, and are used purely for accounting. 

Let S be a subset of service requests. Define the coverage of S by a run R, written cover r(S), 
to be the number of requests in S satisfied by R divided by the total number of requests in S. 
Define the coverage of S by a set U of runs, written coveru(S), to be the average of cover r(S) for 
every run R € U. 

Proposition 3.1 (Average Coverage) Let {S±, S2, S3, ...S a } be a collection of sets of service 
requests such that [J i Si gives all the requests serviced by R* on untrimmed windows. Let U be a 
set of service runs. Ifmhii{coveru(Si)} = fi, then there is at least one service run R £ U such that 
profit(R) > fi-profit(R*). 

The Average Coverage Proposition formalizes the following intuition. Let a group of service 
runs achieve some coverage over a set of requests. Let us also say that we have divided those 
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requests into many different subsets, some of which may overlap, but the union of all the subsets 
is the original set of requests. If we take the subset of requests with the worst average coverage, 
some service run in the group covers a fraction of total requests no smaller than that worst average 
coverage. Otherwise, the average coverage of all subsets would be worse than the coverage of the 
worst covered subset, which is a contradiction. 

Given a way of dividing requests into subsets, we wish to prove that some set of service runs 
achieves a particular lower bound on average coverage. For the case of unit-length time windows, 
we divide requests into subsets based on membership in L, T, or E sets (or smaller partitions 
of these three sets). Since our dynamic programming approach finds an optimal run on trimmed 
windows, to within a factor of 7, we know that this run must perform at least as well as the bound 
we prove. 

For analyzing s = 2 and s = 3, we use ensembles of runs. We average the coverage over sets 
L, T, and E by summing the lowest coverage for the requests in each set and dividing by the total 
number of runs used. By the Average Coverage Proposition, this average gives a lower bound on 
the performance of an optimal run on trimmed intervals. 

Theorem 3.1 Running a repairman algorithm from \12$ at a speedup of s = 2 gives a 27- approximation 
to repairman with unit-time windows in T(n) time. 

Proof: To analyze the performance of the resulting service run, we consider runs A and ^4^, as 
defined above. Refer to Figure [1] for an example. Run A covers all of the requests in set L and all 
of the requests in set T which fall in even numbered periods with respect to R*. Run A R covers 
all of the requests in set E and all of the requests in set T which fall in odd numbered periods. 
For each run, we record a 1 in each set for every period which is completed covered. Note that 
L even is that subset of requests in set L that were serviced by R* in even numbered periods and 
L dd is that subset of requests in set L that were serviced by R* in odd numbered periods, and 
similarly for T and E. We sum these contributions into yields and then divide the yields by the 
total number of runs to get the average coverage for designation. Table Q] illustrates this process, 
in which we get an average coverage of 1/2. Let R be the service run produced by the appropriate 
repairman algorithm from [12j . By the Average Coverage Proposition, profit (R) > max{ profit(A), 
profit(A R )} > (l/2)profit(R*)/~/, giving a 27-approximation. □ 





Leven 


L dd 


± even 


T dd 


Eeven 


E dd 


Run A 


1 


1 


1 











Run A R 











1 


1 


1 


Yield 


1 


1 


1 


1 


1 


1 


Average Coverage 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 



Table 1: Averaging of runs A and A R with speedup s = 2. 



Theorem 3.2 Running a repairman algorithm from ]12$ at a speedup of s = 3 finds a 47/3- 
approximation to repairman with unit-time windows in T(n) time. 

Proof: To analyze the performance, we use four runs: A, A , A, and A R , where A and A R are 
shifted versions of A and A R . Run A follows the same pattern as A but starts the pattern at time 
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1/2 at the location R* has at time 0. Run A R follows the same pattern as A R but starts the pattern 
at time 1/2 at the location R* has at time 1. Run A, now running with speedup s = 3, covers all of 
the requests in L and T but only half the requests in E. Run A R covers all of the requests in T and 
E but only half the requests in L. The shifted runs are copies of A and A R , covering all of T and the 
other half of the requests of E and L, respectively. As before, we sum the contributions and divide 
by the total number of runs to find the average coverages shown in Table [2j The worst coverage for 
any set, averaged over runs A, A R , A, and A R , is 3/4, corresponding to sets L and E. Let R be the 
service run produced by the appropriate repairman algorithm from [12]. By the Average Coverage 
Proposition, profit(R) > max{ profit(A), profit(A R ), profit (A), profit(A R )} > (3/4)profit(R*)/j, 
giving a 47/3-approximation. □ 





Leven 


L dd 


- L even 


T dd 


Eeven 


Eodd 


Run A 


1 


1 


1 


1 


1 





Run A 


1 


1 


1 


1 





1 


Run A R 





1 


1 


1 


1 


1 


Run A R 


1 
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4 
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3 


Average Coverage 


3/4 


3/4 


1 


1 


3/4 


3/4 



Table 2: Averaging of runs A, A, A R , and A R with speedup s = 3. 



4 Canonical Collections and Their Use in Algorithms 

If speedup is a rational number, let s = q/r, where q and r are relatively prime. Otherwise, if the 
speedup is an irrational number, we assume r to be infinite. Let m be the number of requests. 
Below we will define algorithm SPEEDUP such that it runs a repairman algorithm on each set of 
periods from among a collection of min{r, m} different sets of periods. If s is rational and r < m, 
we create sets of periods starting at r time instants spread uniformly from time up to time 1/2. 
If any window would start on any period or half period boundary, we shift the periods by a small 
amount to avoid this case. When r > m, SPEEDUP will use a canonical collection of sets of periods 
which is equally as effective as the r sets (whenever r is finite) but has size just linear in m. The 
procedure CANONICAL defined below "normalizes" the beginnings of all the time windows into 
the range (0,1/2) and then picks a time between each consecutive pair of beginnings to be the 
starting point for a new set of periods. 

Lemma 4.1 Any given set of periods is equivalent to one of the at most m sets of periods in the 
canonical collection, and this collection can be identified in 0(m log m) time. 

Proof: Without loss of generality, we assume that all time windows have unique starting times after 
normalization. Label the normalized starting times of the windows hi, /12, /13, ■ ■ ., h m in ascending 
order as produced by the procedure CANONICAL. The first point chosen by CANONICAL is 0. 
The second point chosen is (hi + /i2)/2, the third is (/12 + h^)/2, and so on. 

A starting point uniquely determines a set of periods, with each successive period starting at 
.5 time units after the start of the previous one. Each set of periods induces a trimming of time 
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CANONICAL 



For time window Sj, 

Let t(si) be the starting time of Sj. 

Let a be the largest integer such that a/2 < i(sj). 

Set hi to be t(si) — a/2. 
Sort the hi values, relabeling them hx, hi, ■ ■ ■ , h m in increasing order. 
Let the first set of periods start at time 0. 
For i from 2 to m, 

Let the i th set of periods start at time (/ti-i + hi)/2. 



windows according to the trimming procedure given before. Define an equivalence relation on sets 
of periods: Consider two trimmings induced by two different sets of periods. If the difference in 
starting times from one set of periods to the next is the same for all time windows, then those two 
sets are equivalent. Consider the time window starting at hi. Those sets of periods are equivalent 
that have starting times strictly between hi and ht+x. Only by moving the starting time earlier 
than hi will the trimming of the earlier window (starting at hi) change, and only by moving the 
starting time later than hi+x will the trimming of the later window (starting at hi+x) change. 
Because our formulation of the traveling repairman problem does not specify a starting position or 
time, all equivalent sets of periods are functionally identical. This equivalence relation partitions 
all possible sets of periods into equivalence classes. By choosing a representative from each of these 
m equivalence classes, we find a collection of m sets of periods that covers all possible trimmings. 

The running time for CANONICAL is 0(m log m) because the m starting times are sorted and 
all other work is linear. □ 

With CANONICAL defined, the description of SPEEDUP using speedup s is straightforward. 
We assume that no window starts on a period boundary as this situation can be easily remedied. 
Let REPAIR be a basic repairman algorithm run after the trimming process has taken place, such 
as the algorithms in [12J. As shown in |12j . REPAIR on a trimmed instance has an approximation 
ratio of at most 7 and a running time of T(n). 



SPEEDUP 

If s is rational and s = q/r in reduced form, where r < m, 
For i from 1 to r, 

Let hi be (i — l)/r. 
Otherwise, 

Let hx, hi, h%, ... h m be the period starting times found by CANONICAL. 
(Note that the total number of hi values is min{r, m}.) 
For i from 1 to min{r, m}, 

Start a set of periods at time hi and trim each window into the period it 

completely fills. 

Run REPAIR on the trimmed windows and retain the best result so far. 



Since algorithm SPEEDUP runs REPAIR a total of min{r, m} times, the running time for 
SPEEDUP is 0(min{r,m}r(n)). 
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Observation 4.1 Algorithm SPEEDUP works for any real number 1 < s < 4; however, our 
analysis will assume that s is a rational number such that s = q/r. In the case that s is irrational, 
our analysis still holds because the functions of s with which we bound performance are piecewise 
smooth and continuous. Thus, we can choose rational numbers arbitrarily close to any real number, 
and our analysis will be well-behaved in the limit. 

5 Analysis of Speedup in the Range 2 < s < 4 

In this section, we describe an ensemble of runs and the coverage given by such a set of runs when 
speedups s is in the range 2 < s < 4. Based on the average coverage of the ensemble, the algorithm 
SPEEDUP produces a 47/s-approximation in this speedup range. 

We assume that no time window starts at time i/{2r) for any integer i. If any do, we can 
perturb all times slightly so that this is not the case. Let [w, w + 1) be any time window. Let u 
be the smallest integer multiple of l/(2r) that is greater than w. We designate subintervals [w,to), 
[w,w + l/(2r)), [w + l/(2r),w + 2/(2r)), [w + (2r-l)/(2r),w+l) by w , w 1} w 2 , w 2r . When 
a given time window is trimmed, some portion of the original window will lie in an earlier period, 
half of it in the period it was trimmed into, and the remaining portion in the following period. As 
before, the goal is to use an ensemble of runs to cover every part of each time window so that we 
can show a lower bound on profit with an averaging argument. However, using sets L, T, and E 
alone is too coarse-grained for general s. Instead, we record the coverage of each subinterval Wi. 
For accounting purposes, we assign a 1 for any subinterval which is covered every period, a 1/2 for 
any subinterval covered every other period, and a otherwise. This accounting scheme implicitly 
averages runs A and A R with their respective shifted versions A and A R . Because the union of 
requests serviced during these subintervals consists of all the requests serviced by an optimal path, 
the Average Coverage Proposition still applies. 

We partition sets L, T, and E into subsets: L\ through L r , T\ through T r , and E\ through E r 
respectively. If we partition each period into r equal- length divisions, then subset Lj, Tj, or Ej, 
for any given j, will consist of those service requests from set L, T, or E, respectively, which were 
serviced by the optimal run in the j division of a period. Note that, for a given period starting 
time, the window for a service request can only overlap with 2r + 1 of these divisions. Thus, we 
can make a one-to-one mapping of the 2r + 1 subintervals of a window to 2r + 1 of the 3r subsets 
of L, T, or E. For a given window, some set of periods maps all subintervals of the window to L\ 
through L r , T\ through T r , and E\ . Another set of periods maps all subintervals of the window to 
L 2 through L r , T\ through T r , and E\ through E 2 . Each remaining set of periods maps one fewer 
Lj subset and one additional Ej subset, where j = i + 1 for set i. Table [3] shows which subsets are 
mapped to subintervals of some time window depending on the set of periods used. 

Figure [2] gives three examples of type A runs for various speedups in the range 2 < s < 4. 
Note that all of the runs in this figure arrive at the same position at the beginning of every second 
period, namely at times 0,1,2, and so on. Portions of runs servicing requests in L, T, and E or 
various subsets are identified: subset E\ being those service requests of E serviced by R* in the 
first half of a period for s = 7/2, subsets T%, T 2 , T3, E%, E 2 , and E3 being those service requests 
of T or E serviced by R* in various quarter periods for s = 11/4, and subsets T%, T 2 , E±, and E 2 
being those service requests of T or E serviced by R* in various thirds of periods for s = 7/3. 

Algorithm SPEEDUP runs a repairman algorithm on min{r, m} different sets of periods and 
takes the best run from those choices. 

Theorem 5.1 For any s in the range 2 < s < 4, SPEEDUP finds a 47 / s- approximation to 
repairman with unit-time windows in 0(min{r, m}T(n)) time. 
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Set of 
Periods 


w 




W 2 . 


w r 


Subintervals 

W r +l W r+2 . 


■ W 2r - 2 


W 2r -l 


w 2r 





Li 




L 3 ■ 


■ T x 


T 2 


T 3 . 


■ T r -l 


T r 


Ei 


1 


L 2 


T 3 


L A . 


■ T 2 


T 3 


T 4 . 


. T r 


Ei 


E 2 




r - 2 


L r -i 


L f 


Ti . 


T r -i 


T r 


Ei . 


E r -3 


E r - 2 


E r -i 


r — 1 


L r 


Ti 


T 2 . 


. T r 


Ei 


E 2 


E T - 2 


E r -i 





Table 3: Possible ways to map the subintervals of a time window to subsets. 
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Figure 2: Examples of type A runs for three different speedups in the range 2 < s < 4, namely at 
7/3, 11/4, and 7/2. 

Proof: For any window [w,w + 1), on service run A, if a trimmed interval starts with subinterval 
Wi, where 1 < i < r, then the window's first q — 2r + i subintervals will be covered each period, and 
the next r subintervals after that will be covered every other period (but no further than the last 
subinterval). When s < 3, run A completely covers subsets Li through L r , covers those requests 
in subsets T\ through T r serviced by R* in half of the periods, and covers those requests in subsets 
7i through T q ^ 2r and subsets E\ through E q _ 2r serviced by R* in the other half of the periods. 
When s > 3, run A completely covers subsets Li through L r and subsets T\ through T r , covers 
those requests in subsets E\ through E q _ 3r serviced by R* in half of the periods, and covers those 
requests in subsets E\ through E r serviced by R* in the other half of the periods. A similar analysis 
applies to A, A R , and A R . 

We assign a 1 for any subinterval which is covered every period and a 1/2 for any subinterval 
covered every other period. For runs A and A together, Wi earns a 1 for all r sets of periods where 
< i < q — 2r, giving a total of r for each such i. For each i > q — 2r, the total decreases by 1/2 
from the total for i — 1. For runs A R and A R together, Wi gets 1 for all r sets of periods where 
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2r — (q — 2r) < i < 2r, giving a total of r for each such i. For each i < 2r — {q — 2r), the total 
decreases by 1/2 from the total for i + 1. 

Now, we take the total over all sets of periods for both runs. For runs A and A together we get 
a total of r for w$. For runs A R and A R together we get a total of r — (1/2) (2r — (q — 2r)) = q/2 — r. 
Summing these together, we get a yield of q/2 for wq. By symmetry, the total for w^ r is also q/2. 
In the case that s < 3, the contributions for runs A and A together decrease by no more than 
1/2 as the contributions for runs A R and A R together increase by 1/2 for each increment of i in 
the range < i < r. In the case that s > 3, contributions from A and A together are constant 
and contributions from A R and A R together only increase or stay constant for < i < r. Thus, 
the yield for all other Wi in all cases is at least q/2. Since r sets of periods for the two pairs 
of runs cost a total of 2r sets of periods to average over, the yield for the algorithm is at least 
q/2 ■ l/(2r) = q/(4r) = s/4. Multiplying the reciprocal by 7 gives a 47/s-approximation. By 
Observation 14.11 this result holds for all real values of s. Note that the details of this proof can 
also be generated using the CREATE- TABLE procedure given in Sect. [6j □ 

6 Analysis of Speedup in the Range 1 < s < 2 

A different analysis gives us better performance as a function of s in the range 1 < s < 2. Based on 
the average coverage of the ensemble, the algorithm SPEEDUP produces a 67/ (s-l-l)-approximation 
in this speedup range. We consider four runs (and their shifts) and balance their relative weightings. 
For the range 1 < s < 2, we represent speedup s as s = (r + k)/r with integers r > 1 and < k < r. 

Our analysis will require several versions of A that have different starting points. To simplify 
notation, let a hop be the distance traveled in time l/(2r) at unit speed. Let A& be run A starting 
A hops further along and run A^ be run A R starting A hops further back. Run A& follows the same 
pattern of movement as run A but starts at t = at the location that R* has at t = — 1/2 + A/(2r). 
Its reverse follows the same pattern of movement as A R but starts at t = at the location that 
R* has at t = 1/2- A/(2r). 

Examples of runs A, A R , A r _k, and A R _ k are shown in Figure [3] for s = 5/4 where r = 4 and 
k = 1. Portions of runs servicing requests in L, T, and E or various subsets are identified: subsets 
L4, Ti, T2, T3, T4, and E\ being those service requests of L, T, and E serviced by R* in various 
quarter periods. We also define A, A R , A r _k, and A R _ k as indicated earlier and implicitly average 
the contributions of these shifted runs with the contributions of their unshifted versions in Tables 
E and El 

In order to create these tables, we define a procedure called CREATE- TABLE that describes the 
process of determining coverage for a particular speedup s for a particular run moved A hops. Note 
that CREATE- TABLE is not an algorithm that is run in the process of finding an approximation 
to a repairman problem with speedup. Rather, it provides a template that can be used to produce 
the tables used in analyzing the performance of such approximations. 

Before CREATE-TABLE can be completely defined, it is necessary to explain the pattern of 
coverage generated by a run. For all type A runs, this pattern takes one of two forms. Let s be 
a rational number such that s = q/r. Type A runs repeat every two periods and thus can be 
represented with a pattern of coverage that uses a 1 to signify a subinterval covered every period 
and a 1/2 to signify a subinterval covered every other period. 

Observe the movement of type A runs, noting that such a run moves forward in one period 
the same distance that run R* moves in q subintervals. During the next period of time, a type 
A run moves backward the same distance that run R* moves during q — r subintervals and then 
forward the same distance that run R* moves during r subintervals. Then, the pattern repeats. 
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Figure 3: Examples of runs A, A R , A r _k, and A R _ k with a speedup of s = (r + k)/r = 5/4. 

When s < 2, run A during its first period in the pattern covers q successive subintervals as it 
moves forward, while in its second period covers r subintervals as it moves forward. Note that 
those subintervals covered as A moves backward add nothing additional to the coverage. Thus, this 
pattern of coverage is represented as r repetitions of 1 and q — r repetitions of 1/2. 



CREATE- TABLE(hops A) 



C(i) 



Let the first element of the coverage pattern be indexed at 0. 
Let the subintervals of interest be numbered through 2r. 
Define function C based on the coverage pattern, such that: 
1 if term i of the pattern is 1 
1/2 if term i of the pattern is 1/2 
otherwise 
Define function T on integers i, where < % < 2r: 

^) = E J r =o^+J-A). 
Define function J- R on the same domain: 
= T(2r - i). 

The final coverage function defined by the table is given by J~{i) + J- R { 



CREATE- TABLE defines functions and such that the final coverage that we define 

is J-{i) + J-{i) R for i = 0, 1, . . . , 2r. Note that the functions J- and J- R are piecewise linear func- 
tions with ranges dependent on the fundamental pattern of coverage. Due to its construction, the 
combination J-{i) + J- R {i) is also a piecewise linear function and symmetrical. Thus, only the range 
< i < r need be listed in tables. 

Although CREATE- TABLE gives a procedure for creating a table for a given speedup, we need 
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tables expressed symbolically to prove coverage for a range of speedups. Instead of using specific 
numbers, we can leave the basic patterns of subinterval coverage for a given style of run (such as 
type A runs) with its shifted version in terms of q and r. By shifting this pattern r times and 
summing the results together, we account for the different alignments a time window might have 
with respect to the various subintervals. This sum is the function which we can now express as 
a piecewise linear function parameterized by q and r. Function T R , which describes type A R runs, 
can be similarly described. To combine the two functions symbolically, we sort the end points of 
the subinterval ranges from both descriptions together. If, for the given range of speedups being 
considered, there are two end points which cannot be ordered, we subdivide the range of speeds 
so that, in each new speed range, the two end points in question can be ordered. Once the end 
points of each subinterval range have been sorted, combining the descriptions from the normal and 
reversed functions of the runs is achieved by simply summing each range. 



r-\i 
r+\k- 
\r + \k 



< i <k 
k < i < r 
r < i < r + k 



F R (i) 



2^ 2 T 2^ 



\k 



r + i 



2r 



k < i < r 

r < i <2r — k 

k < i <2r 



Table 4: Separate coverage functions for runs A and A R when 1 < s < 2. 

In Table H] we use A runs to give an example of how symbolic coverage works. We are only 
interested in the range < i < r. In this range, the sub-ranges < i < k and k < i < r for 
T overlap with the sub-ranges < i < r — k and r — k < i < r for T R . When k < r — k, for 
k < i < r - k, + T R (i) = (r + k/2 - i) + = r + k/2 - i, as in Table When k > r - k, 
for r - k < i < k, + T R (i) = (r - i/2) + (i/2 - r/2 + k/2) = r/2 + k/2, as in Table El The 
combined coverages of runs A, A R , A r _j t , A R _ k , and all of their respective shifted versions are all 
given in Table [5] when k <r — k and in Table [6] when k >r — k. 



Combined contributions 
for A and A R 



Combined contributions 
for and A R k 



r-\i 




o < 


i 


< k 


r + \k - 


i 


k < 


i 


< r 


lr + k- 


¥ 


r — k < 


i 


< r 


k + ^i 




< 


i 


< k 


\k + 2i 




k < 


i 


< r 


¥-k + 




r — k < 


i 


< r 



Table 5: Combined contributions of runs A, A R , A r _k, A R _ k , and all of their respective shifted 
versions, for 1 < s < 2 for the ith interval when k < r — k. Does not include values for i > r 
because those values are symmetrical. 

Lemma 6.1 When we weight the combined contributions from A and A R twice as much as the 
combined contributions from A T _k and A R _ k , the yield for all intervals is at least 2r + k. 
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Combined contributions 
for A and A R 



Combined contributions 
for A r _ fc and A R _ k 



r-\i 




U < 


i 




r — 


A: 


i ,ii 

2 r "i" 2 K 




r — k < 


i 


< 


7 

K 




1 i 1 


i • 

5* 


k < 


i 


< 


r 




k + |i 




< 


i 


< 


r — 


A; 


2' 2 ft 




r — k < 


i 


< 


k 




§r- fc + 




k < 


i 


< 


r 





Table 6: Combined contributions of the same runs as in Table [SJ when k > r — k. 

Proof: We consider yields only for subintervals Wi with < i < r because the cases with r < i < 2r 
are symmetrical. We first handle the case when k <r — k, consulting Table 03 

If < % < k, then the yield for io, is 2r + k + i/2, which is at least 2r + k, since i > 0. 

If k < i < r — k, then the yield for Wi is 2r + 3k/2, which is greater than 2r + k. 

If r — k < i < r, then the yield for W{ is 5r/2 + k — i/2, which is at least 2r + k, since i < r. 

We now handle the case when k > r — k, consulting Table EJ The algebra for the cases when 
< % < r — k and k < i < r gives exactly the same results as the first and third ranges from the 
previous part of the proof. If r — k < i < k, then the yield for subinterval i is 5r/2 + k/2, which is 
at least 2r + k, since r > k. □ 

Theorem 6.1 For speedup s in the range 1 < s < 2, algorithm SPEEDUP finds a 67/(5 + 1)- 
approximation to repairman with unit-time windows in 0(min{r, m}T(n)) time. 

Proof: By Lemma 16-H the yield is at least 2r + k. Since we use two copies each of A and 
A R and a single copy each of A r _j t and A R _ k averaged over r different sets of periods, we apply 
the Average Coverage Proposition over 6r runs. Thus, the approximation ratio is no more than 
6rj/(2r + k) = 6vy/{(r + k) + r) = 67/(5 + 1). By Observation 14. 1\ this result holds for all real 
values of s. □ 



7 Conclusion 

As shown in [13], the technique of trimming provides a powerful tool in attacking route-planning 
problems with time windows. With trimming, it is possible to simplify the structure of many such 
problems and apply an ordering that allows dynamic programming to work. For unrooted prob- 
lems, the cost of this additional order is at most a constant reduction in the profit a run can earn. 
In this paper we have shown that this reduction in profit can be offset, in part or in whole, by 
allowing some corresponding increase in speed over a hypothetical optimal benchmark. Although 
we have only given the techniques needed for unit-time windows here, the extension to the more 
complicated case of window lengths over a range of values is forthcoming in [14J. 

Acknowledgments: We thank Nelson Uhan for bringing a related reference on scheduling to our 
attention. 
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